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We report that, unlike most inclusive quarkonium-production processes, the production rate for 
pp — > J/ip + T + X dominantly depends on two not well-known long-distance nonrelativistic QCD 
matrix elements, (Og ( 3 Si)) and {Oj ( 3 Si)} at leading order in the strong coupling constant, which 
account for the transition probabilities of the color-octet spin-triplet heavy-quark-antiquark pairs 
ccg( 3 5'i) — > J/ip and bbg( 3 Si) — > T, respectively. With the integrated luminosity ~ lOOfb -1 at the 
center-of-momentum energy y^s = 14TeV we expect that approximately 1900 pp — > J/ip + T + X 
events can be observed at the CERN Large Hadron Collider by tagging muon pairs, which are 
enough to probe to the color-octet mechanism. The forthcoming measurement may provide a useful 
constraint to resolve the decade-old puzzle for the polarization of prompt J/ip at the Fermilab 
Tevatron. If corresponding measured rate is significantly less than the prediction, it may imply that 
the current values for the color-octet matrix elements are overestimated. 

PACS numbers: 12.38.-t, 13.85.Ni, 14.40.Pq 
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As an effective field theory of QCD, the nonrelativistic 
QCD (NRQCD) factorization approach [1] has achieved 
great progress in understanding the production and de- 
cay mechanism of heavy quarkonia [2]. The factoriza- 
tion has been proved for the electromagnetic and light- 
hadronic decays [1] and also for a few specific exclusive 
production processes very recently [3-5]. The factoriza- 
tion conjecture for the inclusive production has been em- 
ployed to resolve the large surplus of ip(2S) and prompt 
J/ip at the Fermilab Tevatron [6], which lead to a re- 
markable prediction that the prompt J / ip must be trans- 
versely polarized at the large transverse momentum (jjt) 
[7]. However, the empirical data for the polarization of 
prompt J/ip and ip(2S) measured by the CDF Collabo- 
ration [8, 9] are in disagreement with these predictions. 
That is what we call the puzzle for the polarization of 
prompt J/ip at the Tevatron. 

These predictions strongly depend on the determi- 
nation of the long-distance NRQCD matrix element 
(Of ( 2s+1 Lj)), which accounts for the probability of the 
heavy-quark-antiquark pair QQ with the spectroscopic 
state 2s+1 Lj to evolve into a heavy quarkonium H. 
The color-singlet NRQCD matrix elements (Of ( 3 Si)) 
for H = J/ip and T are well known through pre- 
cisely measured values for their leptonic decay rates. 
On the other hand, the color-octet counterparts such 
as (Of ( 3 5x)), (Of (%)), and (Of ( 3 P 7 )) for J = 0, 
1, 2 are not determined accurately, because in most 
of the hadroproduction rates for H all of these color- 
octet matrix elements involve simultaneously and the 
Pt spectra are not well distinguishable to be fit to the 
data [10, 11]. To make the situation worse, large next- 
to-leading-order corrections dramatically modify the pr 
spectra of these processes [12-15]. As a result, the matrix 
elements determined from various processes are not con- 



sistent with each other [15-21]. In particular, the actual 
values for the color-octet matrix elements (Og^( 1 S'o)) 
and (Og // *( 3 Po)) might be significantly smaller than the 
naive estimates based on the velocity-scaling rules of 
NRQCD [21]. In this situation, it would be highly de- 
sirable to find a hadroproduction process that involves 
fewer numbers of color-octet channels. 

In this paper, we propose that pp — >• J/ip + T + X 
at the CERN Large Hadron Collider (LHC) could be a 
clean probe to the color-octet mechanism of NRQCD and 
compute its pt spectrum and the total production rate 
at the center-of-momentum energy t/s = 14TeV. Unlike 
most inclusive quarkonium-production processes [7, 22- 
24], the production rate for pp — > J/ip+T+X dominantly 
depends on two not well-known long-distance factors 
(0 8 7/V '( 3 5'i)) and (Oj( 3 S±)) and the color-singlet contri- 
bution is suppressed at leading order (LO) in the strong 
coupling constant a s in the asymptotic limit pr — > oo. 
The forthcoming measurement may provide a useful con- 
straint to resolve the decade-old puzzle for the polariza- 
tion of prompt J/ip at the Tevatron. 

This paper is organized as follows. We first summarize 
the NRQCD factorization formula for pp — > J/ijj + T + X 
and describe a systematic way to determine the domi- 
nant contributions. Next we provide detailed arguments 
to convince the dominance of (0%^ ( 3 S\)} and (Oj ( 3 5i)) 
contributions in the px spectrum at LO in a s and the typ- 
ical velocity vq of the heavy quark Q in the quarkonium 
rest frame. Our prediction for the production rate and 
Pt spectrum at the LHC follows and then we conclude. 

The NRQCD factorization formula for the differential 
cross section da of pp — > J/ip + T + X has the following 
schematic form: 

da = f i/p ® f J/p (8 da l3 ^ c - Cni+bln2 (0^)<0^), (1) 
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FIG. 1. Typical Feynman diagrams for gluon-initiated parton 
processes for pp — > J/ip + T + X at order a*. 
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FIG. 2. Typical Feynman diagrams for the color-singlet con- 
tribution to pp — > J/ip + T + X at order a®, which is leading. 



where fi/ p is the parton distribution function (PDF), da- 
is the parton-level cross section, which is perturbatively 
calculable in powers of a s , (O^.) is the NRQCD matrix 
element for the quarkonium Hi involving the NRQCD 
four-quark operator [1]. The subscript n, repre- 

sents the spectroscopic state of the QQ pair. For exam- 
ple, QQm — QQn( 2s+1 Lj) stands for the color-singlct 
(n = 1) or -octet (n = 8) QQ pair with the spin s, 
the orbital angular momentum L, and the total angu- 
lar momentum J for Q = c or b. The symbol ® in- 
dicates the convolution over the partons' longitudinal 
momentum fractions and the summation over partons 
i and j is assumed. According to the velocity-scaling 
rules of NRQCD [1], the production rates for the spin- 
triplet S'-wave quarkonium through the color-octet states 
<2<2 8 ( 3 Si), QQafSo), and QQ 8 { 3 Pj) for J = 0, 1, 2 are 
suppressed by Vq, Vq, and Vq, respectively, compared to 
the color-singlet state QQi( 3 Si). Let us call the product 
of these factors for a process the "velocity-scaling factor" 
V. Therefore, Eq. (1) is a power series in a s , v c , and Vb- 

At LO in ex. g, only gg fusion and qq annihilation con- 
tribute to the parton processes ij — > cc ni + bb n2 of or- 
der cx A s and the channels with = gq and gq are 
missing. At the center-of-momcntum energy of the LHC 
one probes the small- a; region of the PDF, where the 
gluon contribution dominates over the quark contents. 
Therefore, we consider only the gg initial states. Typical 
Feynman diagrams for the gluon-initiated parton process 
gg — > cc ni + bb n2 are shown in Fig. 1, where at least one 
QQ pair is in a color-octet state. Leading color-singlet 
contribution cci( 3 Si) + bbi( 3 Si) appears only at order 
o.g, which are shown in Fig. 2. 

We first consider the color-octet contribution 
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cc s { 3 Si) + 66 8 ( 3 5'i) in Figs. 1 (a)-(d) with V 
which 36 Feynman diagrams contribute. At large px the 
double-fragmentation contribution [Fig. 1 (a)] dominates 
because of the kinematic enhancement. Figure 1 (b) 
also represents mixed contributions ccg( 3 5i) + bbi( 3 Si) 
with V = v 4 and cci( 3 Si) + bb 8 { 3 Si) with V = v 4 , each 



of which has 6 Feynman diagrams. If pt is not large 
enough, then these mixed contributions must dominate 
over cc$( 3 Si) + bb%( 3 Si) by the enhancement factors 1/v 4 
or 1/v 4 while the double- fragmentation contribution 
dominates at large px- 

In addition, cc 8 ( 3 S*i) + bb 8 ( 2s+1 Lj) and cc 8 ( 2s+1 Lj) + 
bb s { 3 Si), where 2s+1 Lj = 1 S or 3 Pj, contribute to dia- 
grams in Figs. 1 (b)-l (d), with V = v" . Here, a, f3 = 3 
for 2s+1 L,j = 1 S and 4 for 3 Pj. The V is comparable 
to that of the double-fragmentation, while they are sup- 
pressed at least by either w 3 or vf in comparison with 
c<5i( 3 Si) + bbi( 3 Si). Although the single- fragmentation 
channel in Fig. 1 (c) may grow up at large pr, that 
contribution is dominated by the double-fragmentation 
[Fig. 1 (a)] by a factor of (m c /pT) 4 or (mb/pr) 4 - Hence, 
it is consistent to ignore ccs{ 3 Si) + bb%( 2s+1 Lj) and 
cc s ( 2s+1 L,j) + bb s ( 3 S!) with 2s+1 Lj = 1 S or 3 P ; over 
the whole pr range. 

As the last color-octet contribution, there are 
cc s { 2s+1 L,j) + bh{ 2s '+ 1 L[ JI ), where 2s+1 Lj, 2s ' +1 L'j, = 
1 So or 3 Pj, whose typical Feynman diagram is shown 
in Fig. 1 (d). In this case V = v"v^ with a, (3 = 3 for 
2s+1 L,j = 1 S'o and 4 for 3 Pj. Thus the contributions are 
suppressed compared to the color-singlct one. Because 
they do not have any fragmentation contributions, they 
are dominated by the double-fragmentation [Fig. 1 (a)] 
by {m c /pT) 4 {mb/pT) A in the large-p^ region. 

So far, we have shown that, at LO in a s and vq, it is 
sufficient to consider only cc$( 3 Si) + bb$( 3 Si), cc$( 3 Si) + 
bbi( 3 Si), and cci ( 3 Si ) + bb$ { 3 S\) to describe the pr spec- 
trum of pp — > J/ip + T + X. The parton-level differential 
cross sections for these processes are available in Rcf. [25] . 

Based on the formalism that we have described earlier, 
we are ready to carry out the numerical calculation of 
the production rates for J/ip + T + X at the LHC. The 
formula (1) involves the factorization of the PDF and 
the parton-level cross section da with the factorization 
scale For the scale ll, we take the transverse mass 
H = rriT = (4toq + Pr) 1 ^ 2 with the heavy-quark masses 
m c = 1.5 GcV and = 4.7 GcV. At LO in a s on which 
we are working, there are no additional hard jets and, 
therefore, pr and tot are defined unambiguously. We 
employ the CTEQ6L parametrization [26] for the PDF. 
We evaluate cx s by setting the renormalization scale to 
be tot and use the next-to-leading-order formula for the 
running coupling constant to be consistent with Ref. [26] . 

In order to evaluate the production rate (1), wc 
have to know the values for the NRQCD matrix ele- 
ments (0^( 3 Si)) for H = J/ip and T, where n — 
1 or 8. The color-singlet matrix element (Of r ( 3 5'i)) 
is usually determined from the leptonic decay rate of 
H, which is the most precisely measured value involv- 
ing H. Wc quote (O^ ( 3 Sx)) = 1.32 GcV 3 [27] and 
(Of ( 3 5i)) = 9.21 GeV 3 [28]. The color-octet matrix cle- 
ment (O'g^ ( 3 S\)} has been fit to the pr spectrum of the 



3 




FIG. 3. The differential cross section da/dpr for pp — > 
J/ip + T + X at y/s =14TeV in units of pb/GeV as a 
function of pr integrated over the rapidity range \y\ < 2.4. 
The solid, dashed, dashed-dotted, and dotted curves repre- 
sent the total, cci( 3 Si) + bb & ( 3 Si), ce&( s Si) + bh( 3 Si), and 
ccs{ 3 Si) + bbg( 3 Si) contributions, respectively. 

inclusive prompt J/tp production rate at the Tevatron in 
the large-pT region. The matrix element (Oj( 3 Si)) has 
been fit [11] to the Tevatron data and used for the po- 
larization analysis [29]. Various determinations of these 
matrix elements can be found, for example, in Ref. [30]. 
In this work, we use (O^^Si)) = 3.9 x 1CT 3 GcV 3 [7], 
(Oj( 3 Si)) = 1.5 x lfr 1 GcV 3 [30]. 

The prediction for da/dpr of pp — > J/ip + T + X at 
y/s =14TeV integrated over the rapidity range \y\ < 2.4 
is shown as the solid curve in Fig. 3. da/dpx vanishes at 
Pt = and increases rapidly until it reaches the maxi- 
mum value da / dpxlmax = 2.9pb/GeV at pt — 1.5 GeV. 
Then it monotonically decays as pt increases. Near 
Pt = 0, cci( 3 S , i)+668( 3 S'i) dominates and it is the largest 
for pt < 4 GeV. At large pr, the double- fragmentation 
contribution cc 8 ( 3 5i) + 66g( 3 5'i) [Fig. 1(a)] dominates 
and it is the largest for pt > 6 GeV. In the remaining re- 
gion, 4 GeV < pt < 6 GeV, all of the three contributions 
compete together. 

As is mentioned earlier, we have ignored the color- 
singlet contribution in Fig. 2 because it is suppressed 
compared to the color-octet one by a relative order of 
Qj. In order to assure that this approximation is safe, we 
make a rough estimate of the color-singlet contribution 
to pp — > J/tp + T + X. The corresponding parton pro- 
cesses consist of 3 types: Types A, B, and C are shown 
in Fig. 2 (a)-2 (c), respectively. 

We first classify the scaling behavior of the amplitude 
for the double-fragmentation diagram shown in Fig. 1 (a). 
The propagator for the exchanged gluon has the scaling 
1/Pt and the gluon propagators attached to the cc and 
bb pairs are of order 1/m 2 . and 1/m 2 , respectively. The 
product of the two triple-gluon vertices must scale as 



the typical momentum squared p T . Therefore, the resul- 
tant relative scaling of the double-fragmentation process 
is v\v\ I (m^mf) , where we have included the velocity- 
scaling factor V = v\v\. 

In the case of type A, there are four heavy-quark prop- 
agators and one gluon propagator which have a large 
momentum transfer that accounts for the scaling factor 
l/p T - The phase space for the cci + bb\ + gg final state 
enhances the scaling factor by an order of p T in compar- 
ison with that for the two-body final state cc$ + bbg of 
the double- fragmentation process. As a result, the type 
A process is scaled by l/p T compared to the double- 
fragmentation process. In a similar manner, the scaling 
factor for type B is the same as that of type A. In types 
A and B, there are extra hard jets in the final states. 
Therefore, the color-octet contribution, which does not 
have hard jets, is distinguished from these color-singlet 
contributions that can simply be removed by imposing 
an appropriate veto. Type C diagrams involve finite box 
diagrams with two virtual gluons whose momentum must 
be of order the typical momentum transfer pt- A rough 
estimate of the scaling can be found by substituting a 
typical momentum transfer pt to the two gluon propa- 
gators, four heavy-quark propagators, and the measure 
of the loop momentum. The relative scaling of type C 
process is again l/p T . 

Therefore, the suppression factor for the 
color-singlet contribution relative to the 
color-octet double-fragmentation process is 
[a s /(47r)] 2 [m c /(u c p T )] 4 x[m h /(w b p T )] 4 where we have 
included the strong coupling suppression factor a 2 to 
the color-singlet contribution. However, if pt is small, 
then the factors pt in this scaling should be of order tuq 
and, therefore, the suppression factor becomes roughly 
a 2 s / [(Air) 2 vl] in comparison with cc\{ z Si) + bb 8 { 3 Si), 
which dominates over the double-fragmentation con- 
tribution as pt — > 0. This factor is of order 10~ 2 if 
we assume v 2 ~ 0.1 and we use the renormalization 
scale of order mi,. According to this rough estimate, we 
expect that our prediction may not be contaminated 
by an order-a^ color-singlet contribution even at small 
Pt- If we introduce a lower pt cut pt > 5 GeV, where 
the double-fragmentation rises up, then the color-octet 
contribution becomes more significant. In addition, by 
imposing a veto that the final state must not include 
hard jets, one can remove the color-singlet contribution 
of types A and B, which are potential sources of a large 
background. 

The kinematic enhancement of the double gluon frag- 
mentation can also be applied to the color-singlet pro- 
cesses of order a® through the subprocesses gg, qq — » gg 
followed by the fragmentations of the two final-state glu- 
ons into cci( 3 Si) + gg and bbi( 3 Si) + gg, respectively. 
This process is suppressed compared to the color-octet 
double gluon fragmentation g ccg( 3 5i) and g — > 
&6 8 ( 3 S'i) that we have calculated in this paper by a factor 
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of order a^/[(47r) 4 t^t^] ~ 7 x 10~ 5 with an assumption 
a s = 0.2 at threshold with the scale fi = Ivtiq} In addi- 
tion, one can consider the mixed double gluon fragmen- 
tation of order or s through the subprocesses gg, qq — > gg 
followed by the fragmentations of the two final-state 
gluons into ccs(bb 8 )( 3 Si) and bbi(cci)( 3 Si) + gg, re- 
spectively. These mixed double-fragmentation contri- 
butions are also suppressed relative to the color-octet 
counterpart by factors of a 2 / [(At:) 2 v^] ~ 3x 10~ 2 and 
a 2 /[(A-K) 2 v^] — 3 x 1CT 3 for the cc 8 + bbi and bb 8 + cci 
fragmentations, respectively. 2 In summary, the mixed 
double-fragmentation contributions may occupy only a 
few percents of the color-octet double fragmentation and 
the pure color-singlet double-fragmentation contribution 
seems to be well suppressed compared to the color-octet 
double- fragmentation contribution. 

Based on the NRQCD factorization formalism, we have 
computed the px spectrum da/dpr and the total cross 
section atot for pp — > J/tp + T + X. The short-distance 
coefficients were computed at LO in a s and vq. atot in- 
tegrated over the rapidity and the transverse-momentum 
range \y\ < 2.4 and pr < 30GeV is predicted to be 
0tot = 13 pb at = 14TeV. Because the production 
rate is subject to the numerical values for (O'^ ( 3 Si)) 
and (Oj ( 3 Si)} , the forthcoming empirical rate may pro- 
vide information whether the values were reasonable or 
not. 

The feature of this LO process that the color-singlct 
channel is suppressed provides us with a clean probe to 
the color-octet mechanism. In comparison with the in- 
clusive single-quarkonium production at the hadron col- 
liders, the suppression of the color-singlet channel is more 
significant by a relative factor of [a s j (47r)] 2 [mQ / (vqpt)} 4 
at large pr- The pr spectrum is dominated by cci ( 3 5*i) + 
bbs{ 3 Si) and ccs ( 3 Si ) + bbs( 3 Si) channels at the low- and 
large-pr region, respectively, and the cc$( 3 Si) + bbi( 3 Si) 
contribution is comparable to those two listed above only 
around the region 4 GeV < pr < 6 GeV. Assuming the 
integrated luminosity ~ 100fb _1 at <Js =14TeV and 
considering the branching fractions B[J/ip — > fi + /i~] = 
5.93% and B[T -> = 2.48% [32], we expect that 



1 Our order-of-magnitude estimate relies somewhat on the extra 
factors of (4-7r), which arc from the phase space factor or loop 
integration. It is well known that the inclusion of the factor (47r) 
improves the naive estimates, for example, in the muon decay 
rate using the Fermi's four-fermion interaction Lagrangian for 
charged weak interaction. Our naive estimates may break down if 
there are large enhancement factors in the color-singlet channels. 
An explicit calculation of the complete color-singlct contribution 
may reveal if our rough power counting is indeed reasonable. 

2 If we make use of the fragmentation probability at threshold 
given in Ref. [31], we can make rough estimates on the ratios 
of the double-fragmentation processes cci + bbi , ccg + bbi , and 
fobs + cci compared to the ccs + bb~8 fragmentation as ~ 2 X 10~ 4 , 
~ 5 X 10~ 3 , and ~ 4 X 10~ 2 , respectively. These values are 
consistent with the scaling factors listed above. 



approximately 1900, 520, and 160 events can be observed 
by tagging muon pairs under the cuts pr > 0, 5, and 
10 GeV, respectively, at the LHC. If one can improve the 
acceptances for J/ip and T by extensive Monte Carlo 
studies of final-muon pairs, then the observation of the 
events can be quite promising in the near future. Inclu- 
sion of e + e~ decay modes of J/ip and T may increase the 
number of events by a factor of 4. One can also improve 
the prediction by including the subprocess qq — » J/ijj + T 
via two-gluon exchange that is neglected in this work. If 
the forthcoming measured rate is significantly less than 
our prediction, then it may be an indication that the 
current values for the octet matrix elements are overes- 
timated. 

We consider possible sources that may contaminate 
our predictions. In this work, we have computed the 
rate for the direct J/ip's rather than prompt J/^'s that 
include feeddowns from higher resonances and that ex- 
clude those from the B-meson decay. While nonprompt 
signals from B decays can be eliminated with the aid of 
silicon vertex detectors, the separation of the feeddowns 
from higher resonances is not easy. Therefore, the simul- 
taneous determination of the color-octet NRQCD ma- 
trix elements (O'^^Si)) and (Oj( 3 Si)) using this pro- 
cess may be nontrivial. However, the dominance of the 
cc8( 3 <Si) + bbg,( 3 Si) contribution may still be true, even 
after including the feeddowns from higher resonances, 
especially at large px due to the double-fragmentation 
process. Once the px distribution for the J/ip + T pro- 
duction is measured, it may provide us with significant 
constraints on the product of the two matrix elements. 

We note that the total rate and the pr spectrum for the 
inclusive single-quarkonium production at hadron col- 
liders suffer from large next-to-lcading-order corrections 
in a s . In particular, the next-to- leading-order correc- 
tions to the color-singlet contribution to the inclusive 
J/ip production enhance the rate by an order of mag- 
nitude [12-14]. Therefore, it is natural to worry that the 
inclusion of order-a® contributions to this process may 
modify the px spectrum significantly. However, in the 
case of J/tp + T production, the suppression factor of 
the color-singlet channel to the color-octet one is a 2 s /p T 
which is significantly smaller than the corresponding fac- 
tor l/p T for the single-quarkonium production, in which 
the color-singlct and color-octet contributions are of the 
same order (a 3 ) at LO. We anticipate that our argu- 
ments can be tested quantitatively by measurements and 
by explicit calculations of complete order-a^ contribu- 
tions to the color-singlet channel in near future. Once 
the suppression of the color-singlet channel is confirmed 
at that order, then the measurement of the px spec- 
trum of J/ip + T at the LHC may play an important 
role to determine the color-octet NRQCD matrix ele- 
ments (Og^^Si)} and (Oj^S 1 !)). Even in the worst 
case, that one could not observe the pr spectrum at the 
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desired level, at least one can obtain strong upper bounds 
to both NRQCD matrix elements, that were not easy to 
extract empirically in other experiments. Therefore, the 
forthcoming measurement may provide useful informa- 
tion that can be used to understand the mechanism of 
the decade-old puzzle for the polarization of prompt J/ip 
at the Tevatron. 

In summary, we have computed the production rate 
and pt distribution for pp — > J ftp + T + X at the 
LHC in the framework of the NRQCD factorization 
formalism at leading order in a s . Unlike most inclu- 
sive single-quarkonium-production processes, in which at 
least three independent color-octet NRQCD matrix ele- 
ments (Of ( 3 5!)), (Of CSq)), and (Of ( 3 Pj)) involve, 
this process depends dominantly on two matrix elements 
(O^C^i)) and (Oj( 3 Si)) at LO in a„, especially at 
large pr- If the corresponding measured rate is signifi- 
cantly less than the prediction presented in this work, it 
may imply that current values for the octet matrix ele- 
ments are overestimated. The forthcoming measurement 
of pr distribution may probe to the color-octet mecha- 
nism in quarkonium production, by imposing strong con- 
straints to these NRQCD matrix elements. 
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